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MODEL  VERSUS  PATH  DATA 

H.  J.  Hebe,  K.  C.  Allen,  G.  R.  Hand,  R.  H.  Espeland,  and  E.  J.  Violette* 

A  practical  atmospheric  millimeter-wave  propagation  model  (MPM) 
is  updated  and  tested  with  experimental  data  from  horizontal,  line- 
of-sight  links  when  there  is  no  precipitation.  The  MPM  computer 
program  predicts  attenuation  and  delay  properties  of  moist  air  over 
ranges  in  frequency  from  1  to  1000  GHz  and  in  height  from  0  to  30  km. 

Input  variables  are  radio  path  distributions  of  pressure,  temperature, 
relative  humidity,  and  a  suspended  droplet  concentration  simulating 
haze  and  fog  conditions.  Terrestrial  path  data  from  millimeter-wave 
propagation  experiments,  including  those  from  a  27  km  link  operated 
at  11.4,  28.8,  and  96.1  GHz  by  the  Institute  for  Telecommunication 
Sciences  (ITS),  have  been  analyzed.  Calibrated  mean  signal  levels 
permitted  studies  of  water  vapor  losses.  In  addition,  a  spectra) 
analysis  was  performed  of  clear-air  scintillations  caused  by  turbu¬ 
lence.  In  general,  good  agreement  is  obtained  with  the  MPM  for  test 
frequencies  up  to  430  GHz. 

Key  words:  atmospheric  attenuation  and  delay;  millimeter  wave  properties 

of  moist  air;  propagation  program  MPM;  terrestrial  radio  path  data 

1.  INTRODUCTION 

Atmospheric  propagation  limitations  dominate  most  considerations  in  the  advance¬ 
ment  of  millimeter-wave  applications  (Crane,  1981;  Bohlander  et  al.,  1985).  Adverse 
weather  causes  millimeter-wave  signal  degradations  due  to  rain,  wet  snow,  suspended 
particles,  and  water  vapor.  A  propagation  model  provides  a  cost-effective  means  of 
predicting  the  performance  of  a  system  for  its  intended  use  by  considering  limiting 
factors  of  the  atmosphere  that  in  the  actual  operating  environment  may  be  difficult 
to  identify.  The  Institute  for  Telecommunication  Sciences  (ITS)  has  developed  a 
modular  Millimeter- wave  Propagation  Model  labeled  MPM  (Allen  and  Llebe,  1983;  llebe, 
1983).  Based  on  meteorological  variables,  the  model  predicts  propagation  effects 
(that  is,  radio  path  attenuation,  delay,  and  medium  noise)  for  a  radio  path  through 
the  first  100  km  (nominally  ground-to-space)  of  the  atmosphere.  Frequencies  up  to 
1000  GHz  (1  THz)  are  considered.  The  model  has  been  used,  for  example,  to  provide 
estimates  of  millimeter-wave,  average-year  attenuation  distributions  for  diverse 
climates  across  the  United  States  (Allen  et  at.,  1983). 

This  report  addresses  the  moist  air  portion  of  the  MPM.  Central  to  the  model 
are  more  than  450  spectroscopic  parameters  describing  local  (<1  THz)  absorption 
lines  of  the  molecules  0^  and  H^O  complemented  by  continuum  spectra  for  dry  air, 

*The  authors  are  with  the  Institute  for  Telecommunication  Sciences,  National 
Tel  ecou-Mur.icatiors  and  information  Administration,  U.  S.  Department  of  Commerce, 

325  Broadway,  Boulder,  Colorado  80303. 
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water  vapor,  and  hydrosols  (i.e. ,  suspended  water  droplets  for  haze  and  fog  condi¬ 
tions).  Laboratory  measurements  of  absolute  attenuation  rates  at  138  GHz  for  simu¬ 
lated  air  with  water  vapor  pressures  up  to  saturation  led  to  an  improved  water  vapor 
continuum  (Liebe,  1984).  In  addition,  other  research  findings  have  been  reported 
recently  on  the  spectroscopic  data  base,  which  impact  predictions  and  thus  warranted 
an  update  (Liebe,  1985)  of  earlier  versions  (Waters,  1976;  Liebe,  1981).  The  new 
HPH  describes  millimeter-wave  characteristics  of  moist  air  and  is  presented  in 
Section  2.1  in  the  form  of  a  comprehensive  complex  refractivity  N. 

Terrestrial  propagation  experiments  ir.  the  10  to  100  GHz  range  have  been  con¬ 
ducted  by  ITS  over  path  lengths  between  0.1  and  27  km  employing  cw  and  broadband 
(>1  GHz)  signals  (Violette  et  at.,  1983;  Espeland  et  al.,  1984).  In  order  to  com¬ 
pare  MPM  predictions  with  realistic  data,  absolute  intensity  measurements  were  per¬ 
formed  of  phase-coherent  signals  at  11.4,  28.8,  and  96.1  GHz  over  a  27.2  km  path 
located  in  Boulder,  Colorado.  Calibrated  levels  of  received  signal  amplitudes  have 


been  processed  for  time  periods  free  of  precipitation  to  obtain  mean  values  s(dB) 


and  variances  o$(dB  )  as  indicators  for  atmospheric  losses  and  fluctuations  intro¬ 


duced  by  water  vapor,  turbulence,  and  possible  layering  effects.  Attenuation  (dB/km) 
measurements  have  been  made  at  28.8  and  96.1  GHz  when  the  atmosphere  wes  well  mixed. 

Natural  atmospheric  fluctuations  in  density  and  composition  cause  variations  in 
the  refractivity  N±AN  on  various  time  and  space  scales  outlined  in  Table  1. 
Turbulence- Induced  variability  AN  serves  to  explain  temporal  and  spatial  fluctu¬ 
ations  in  received  signal  amplitudes,  phases,  and  angles-of-arrival.  Turbulence 
models  addressing  amplitude  scintillations  of  millimeter  waves  art  discussed  briefly 
in  Section  2.2. 


Table  1.  Characteristic  Time  and  Horizontal  Scales  of 
Atmospheric  Processes  (Orlanski,  1975) 


Process 


Approximate  Time  Scale  Approximate  '.’orizontal  Scale 


Turbulence 

10'2-104 

(millisec.  to 
hours) 

10*3-103 

(micro) 

Earth  rotation 

105 

(day) 

103-104 

(meso) 

Weather  systems 

105-106 

(week) 

104-105 

(synoptic) 

Seasonal  changes 

106-107 

(month) 

105-106 

(macro) 

Weather  cycles 

4-* 

o 

'■N, 

1 

o 

(year) 

106-107 

(global ) 

v-  v  \. 


Experimental  data  capable  of  substantiating  these  models  have  been  limited  to  a  few 
cases  (e.g. ,  Cole  et  al.,  1978;  Filho  et  a!.,  1983;  Vogel  et  al.,  1984).  Scintilla¬ 
tion  results  obtained  with  the  Boulder  link  for  a  frequency  ratio  of  1:8.4,  add  suit¬ 
able  test  data.  Section  2.3  touches  upon  the  problem  that  weak  turbulence  might 
sustain  large-scale  refractive  index  stratifications.  Such  conditions  might  result 
in  deep  fade  events  since  stable  layering  can  lead  to  multiple  rays  that  mutually 
interfere. 

Details  of  the  ITS  terrestrial  propagation  experiment  are  summarized  in  Sec¬ 
tion  3.  Observations  of  received  signal  strengths  obtained  with  this  link,  together 
with  simultaneously  measured  meteorological  data,  are  presented  in  Section  4.  Also 
included  are  selected  data  up  to  430  GHz  that  have  been  reported  in  the  literature. 
Section  5  contains  the  results  of  various  intercomparisons  between  tractable  propa¬ 
gation  data  from  individual  events  and  model  results.  The  report  concludes  w'th  a 
discussion  of  the  validity  of  predictions  and  their  value  in  forecasting  the  perform¬ 
ance  of  millimeter  wave  systems  operating  through  the  atmosphere. 


2.  MILLIMETER-WAVE  CHARACTERISTICS  OF  MOIS1  AIR 
Millimeter  waves  traversing  an  atmospheric  line-of-sight  (LOS)  path  lose  energy 
in  various  ways,  most  of  which  are  frequency-dependent.  Radiated  energy  can  be  ab¬ 
sorbed,  scattered  out  of  the  path,  or  diverted  by  refraction  and  diffraction  effects 
that  cause  the  signal  either  to  miss  the  receiving  antenna  and/or  to  arrive  out  of 
phase  with  the  main  ray.  All  these  interactions  are  dynamic  in  nature,  driven  by 
the  variable  atmosphere  and  leading  to  modulations  of  the  propagating  radio  wave  in 
amplitude  and  phase  (transit  time). 

Two  basic  approaches  for  predicting  transmission  characteristics  of  moist  air 
are  taken: 

’  Physical  Model  (Section  2.1).  Number  concentrations  of  absorbers,  refractive 
layers,  and  diffraction  boundaries  as  part  of  the  path  geometry,  all  are 
accounted  for  based  on  physical  quantities  measurable,  at  least  in  principle, 
at  any  given  path  location  and  time. 

The  transmission  of  broadband  (>100  MHz)  signals  adds  the  additional  dimension 
of  frequency  response  across  the  signal  bandwidth.  In  digital  schemes,  disper¬ 
sive  group  delay  within  the  atmospheric  propagation  channel  can  cause  inter¬ 
symbol  interference  limiting  the  maximum  data  transmission  rate.  For  example, 
100  megabit  per  sec  (Mb/s)  tolerate  about  3  ns  deviations  from  frequency- 1  inear 
behavior  across  the  bandwidth  with  no  increase  in  bit  error  rate  (BEK). 


•  Statistical  Models  (Sections  2.2  and  2.3).  Time  series  of  received  signal  level 
(RSL)  are  analyzed  in  conjunction  with  meteorological  data  over  time  intervals 
of  events  specified  in  Table  1.  The  RSL  is  treated  as  a  random  variable,  dis¬ 
regarding  the  physical  origin  of  the  variability.  Typically,  a  median  value  is 
defined  over  a  reference  period  (e.g. ,  1  h)  superimposed  by  short-term  fluc¬ 
tuations  and  both  processes  are  evaluated  statistically. 

2.1  The  Propagation  Model  MPM 

The  computer  program  MPM  is  a  physical  model  developed  for  applications  in  tele¬ 
communications,  radio  astronomy,  remote  sensing,  etc.,  within  the  1  to  1000  GHz  range 
(Liebe,  1981;  1985).  Radio  path  behavior  through  a  moist  atmosphere  is  described  by 
cumulative  path  attenuation 

d 

A  =  J"a(x)  dx  dB  (la) 

o 

and  by  path  delay 
d 

B  =  Jp(x)  dx  ps,  (lb) 

o 

where  dx  is  an  increment  of  the  path  length  x  in  kilometers  (km)  which  is  of  total 
length  d.  Attenuation  A  quantifies  the  amount  of  energy  extracted  from  a  plane  wave 
propagating  through  the  atmosphere,  delay  p  is  a  measure  of  the  excess  traveling 
time  with  reference  to  vacuum.  Specific  power  attenuation  a  and  propagation  delay  p 
can  be  expressed  in  the  form 


where  frequency  f  is  in  gigahertz  (GHz)  throughout.  The  heart  of  the  model  is  a 
macroscopic  measure  of  interactions  between  radiation  and  the  propagation  medium, 
expressed  as  complex  refractivity  in  N  units  (10  ^  h  ppm) 


N  -  N  ♦  N'(f)  ♦  JN"(f ). 


(3) 
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The  refractivity  consists  of  a  frequency- independent  tern  NQ  plus  various  spectra  of 
refractive  dispersion  H'(f)  and  absorption  N"(f). 

Meteorological  Parameters 

Gaseous  oxygen  (O2)  and  water  vapor  (i^O),  and  suspended  water  droplets  are 
considered  to  be  the  principal  absorbers.  The  physical  state  of  air  is  described  by 
four  measurables  P-T-RH-w,  which  relate  to  internal  model  variables  p-0-e(v)-wQ  as 
follows: 


P  =  p  +  e 


RH  =  (e/es)100 


(kPa),  T  =  (300/0)  -  273.15  (°C), 


w  =  wQ  G(RH) 


(g/»J). 


Barometric  pressure  is  labeled  P,  where  p  is  dry  air  pressure  and  e  is  partial  water 
vapor  pressure,  all  in  units  of  kilopascal  (1  kPa  =  10  mb);  temperature  T  in  degrees 
Celsius  (°C)  is  converted  to  a  relative  inverse  temperature  parameter  u,  which  for  T 
in  absolute  Kelvin  units  equals  300/T(K);  and  relative  humidity 

RH  =  (e/e$)100  =  (v/vs)100  =  ml.  51  e  e'5  10(9,834  0  "  10)  i  1QCX  (4) 

is  given  as  absolute  humidity  ratio  described  interchangeably  by  either  water  vapor 
partial  pressure  e  or  vapor  concentration  v  =  7.217  e  0  (g/m3).  The  subscript  "s" 

denotes  the  saturation  limit  at  RH  =  100%.  A  dry  mass  concentration  of  hygroscopic 
aerosol  is  wq  in  g/m  and  an  approximate  expression  for  the  unitless  mass  growth 
factor  for  RH  <  96%  is  given  by  (Haenel,  1976) 


G  r  100/(100- RH). 


Complex  Refractivitv  M 


Radio  refractivity  of  air  has  been  measured  accurately  at  microwave  frequencies 
(e.g. ,  Boudouris,  1963)  and  reviewed  recently  (Hill  et  ah,  1962)  to  be 

Nq  =  (2.588  p  ♦  2.39  e)0  ♦  41.63  e  O2.  (6) 

Absorption  and  dispersion  spectra  are  obtained  from  n  -  4S  o\  '  s  ir.vs, 
nfa  =  30  water  vapor  lines,  and  various  continuum  spectra  (dry  air),  {water 
vapor),  and  (hydrosol)  via 


.  ■  •  ••  ;,v 


5 


d  =6.67  D2f  10"3 
o 


is  exceeded  (e.g. ,  f  =  96.1  GHz,  D  =  0.3  m,  dQ  £  0.058  km). 


b)  Free  space  loss, 


Lq  =  20  log  (f  d)  +  92.45  dB,  (28) 

accounts  for  the  decrease  in  rf  power  with  distance  d  (km)  from  an  isotropic  trans¬ 
mitter  antenna  (e.g.,  f  =  96.1  GHz,  d  =  27.2  km,  Lq  =  160.8  dB). 
c)  The  receiver  noise  power  level  equals 

Qr  =  10  log  (Af  Te)  -  138.6  dBm,  (29) 

where  is  the  effective  receiver  noise  temperature  (K)  which  includes  mixer  rf 
loss  and  antenna  contributions,  and  the  detection  bandwidth  Af  is  in  MHz  (e.g., 

Te  =  1200  K,  Af  =  1.0  MHz,  Qr  =  -107.8  dBm). 

To  measure  the  highly  variable  atmospheric  loss  A,  the  following  parameter  have 
to  be  considered  for  the  LOS  experiment: 


d . for  the  propagation  path, 

Pt,  Dt,  Lj. ,  f . for  the  transmitter, 

Tr,  D  ,  L  ,  lu,  Af . for  the  receiver, 


Tc,  D  ,  L  ,  lM,  Af . for  the  receiver, 

t  r  r  M 

in  order  to  ensure,  under  all  conditions,  a  system  performance  of  F  >  0  dB. 

3.1  Link  Description  (27-km  Path) 

The  propagation  path  was  located  in  Boulder,  CO,  and  had  a  length  d  =  27.2  km. 
The  first  7- km  segment  starting  at  the  receiver  terminal  is  across  the  city  of 
Boulder  and  the  remainder  is  over  open  fields  and  farms.  This  path  runs  generally 
parallel  to  the  foothills  on  the  west,  approximately  5  to  7  km  from  the  first  moun¬ 
tain  ridge.  Receiver  and  transmitter  terminal  locations  at  h(RX)  =  1.654  and 
h(TX)  =  1.646  km  elevation,  respectively,  and  a  cross  section  of  the  propagation 
path  are  shown  in  Figure  3.  Antenna  beamwidths  BW  at  the  operating  frequencies  are 
indicated  as  well  as  ray  paths  for  a  standard  (k  =  1.34,  effective  Earth  radius)  and 
a  subrefracti ve  (k  =  0.29)  constant  refractivity  gradient  (e.g.,  Stephansen,  1981), 
where  the  latter  condition  would  cause  the  ray  to  just  touch  the  ground  in  the  mid¬ 
dle  of  the  path. 
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(RSL)  s^  2  3  (dB)  qualified  the  field  experiment  for  testing  predictions  made  with 

propagation  models  discussed  in  Section  2.  Examples  of  daily  RSL  recordings,  free 

of  precipitation,  are  presented  and  a  power- spectral  analysis  of  relatively  quiet 

and  relatively  active  fade  periods  is  performed. 

Together  with  the  three  cw  channels,  a  wide-band  channel  was  operated  using  a 

30.3  GHz  carrier  with  a  500  Mb/s  biphase-shift-key  modulator  and  a  1  ns  resolution 

impulse  probe  (Violette  et  al.,  1983).  Bit-error-rates  (BER)  were  measured  down  to 
”8 

about  2  x  10  for  a  1  second  gate-time  and  a  pseudo-random  code.  In  clear  weather, 
the  error  rates  stayed  usually  below  the  detection  threshold  and  are  of  no  concern 
in  this  report. 

Performance  of  an  LOS  system  is  evaluated  by  the  signal-to-noise  ratio  (fade 
margin) 

F  =  Sy  -  Lo  -  Qr  -  A  dB  (24) 

realized  at  the  receiver.  The  quantities  in  (24)  represent  a  system  gain  S  ,  the 
free  space  loss  Lq,  receiver  noise  power  N  ,  and  atmospheric  propagation  loss  A  (la); 
all  but  A  are  examined  in  more  detail  in  the  following, 
a)  A  system  gain  factor 

Sy  =  10  log  Pt  +  Gt  +  Gr  -  Lt  -  Lp  -  L^  dBm  (25) 

is  the  sum  of  various  hardware  contributions,  where  Pt  is  the  average  transmitter 
power  in  milliwatts,  Gt  is  the  directional  gain  of  the  transmitter  antenna  toward 
the  receiver;  Gr  is  the  directional  gain  of  the  receiver  antenna  tov/ard  the  trans¬ 
mitter;  Lt  and  L^  are  rf  losses  of  transmitter  and  receiver  antenna  feeds,  respec¬ 
tively;  and  L^  is  a  conversion  efficiency  (loss)  of  the  receiver  front-end.  Para¬ 
bolic  dish  antennas,  the  type  typically  used  in  millimeter-wave  systems,  display  a 
maximum  directional  gain 

Ga  =  20  log  (fO)  +  10  log  g  +  20.4  dB  (26) 

where  0  is  the  antenna  diameter  in  meters,  and  q^  is  the  efficiency  factor  for  a 
particular  design  (e.g. ,  f  =  96.1  GHz,  D  =  0.30  m,  q  =  0.6,  G  =  47  dB).  Geometric 
dish  imperfections,  illumination  blockage,  and  aperture-to-medium  coupling  loss  de¬ 
crease  the  efficiency  q.  Full  antenna  gain  G  is  realized  when  a  minimum  distance 
(far- field  range) 


atmospheric  conditions  (weak  turbulence)  plus  rapid,  height-dependent  reductions  in 
water  vapor  concentration  (e.g. ,  fog  over  moist  river  valleys  during  calm,  cool 
nights,  and  early  mornings  of  hot,  summer  weather).  Under  such  conditions  a  radio 
ray  can  be  trapped  and  channeled  in  a  duct,  where  the  normal  spherical  spreading 
loss  L  (see  Section  3)  is  no  longer  valid  and  possible  dispersive  propagation  velo¬ 
cities  exist.  The  result  is  a  strong  short-term  RSL  variability  due  to  more  than 
one  component  of  the  radio  field  arriving  at  the  receiver  by  way  of  separate,  slight¬ 
ly  differing  paths.  Occasionally,  when  only  two  rays  interfere  with  each  other, 
both  being  of  similar  magnitude,  the  RSL  can  behave  very  erratically,  with  unpredic¬ 
table  deep  fades  in  the  range  of  20  to  40  dB  and  serious  amplitude  distortions  across 
the  propagation  channel  bandwidth.  Multipath  effects  present  problems  for  the  trans¬ 
mission  of  digital  data.  Delay  spread  causes  the  received  symbols  to  overlap.  The 
resulting  intersymbol  interference  leads  to  decoding  errors.  For  example,  a  delay 
difference  of  3  ns  can  raise  the  BER  at  data  transmission  rates  greater  than  30  Mb/s  . 
A  channel  transmission  model  capable  of  simulating  three-ray  multipath  fades  was 
given  by  Rummler  (1979). 

Multipath  fading  can  also  be  expected  when  reflections  occur  from  terrain  fea¬ 
tures  such  as  hills,  buildings,  trees,  etc.  The  fade  depth  depends  on  the  path  clear¬ 
ance,  surface  curvature,  surface  roughness,  etc.  A  strong  specular  component  will 
produce  regular  (A/2)  interference  patterns  when  the  receiving  antenna  is  moved  in 
space.  Another  phenomenon  is  obstacle  fading  due  to  weather-dependent  refractivity 
changes.  The  effective  earth  radius  (typically  4/3)  changes  so  that  vertical  clear¬ 
ances  become  marginal  from  large-sized  obstacles  (hills,  the  ground,  buildings,  tree 
tops,  etc.)  close  to  the  LOS  path.  These  fades  are  not  particularly  frequency- 
sensi ti ve. 

It  is  possible  to  reduce  the  effect  of  multipath  fading  in  a  number  of  ways; 
the  most  common  are  narrower  beamwidth,  antenna  tilts,  space  and  frequency  diversity 
and  good  path  clearances. 


3.  BOULDER  LOS  EXPERIMENT 

A  terrestrial  1 ine-of-sight  (LOS)  link  designed  to  investigate  atmospheric  milli 
meter-wave  propagation  (Violette  et  al.,  1983;  Espeland  et  al.,  1984)  was  operated 
over  a  27.2-km  path  near  Boulder,  CO,  between  August  9,  1983,  and  June  14,  1984. 

The  link  was  outfitted  at  the  receiver  site  with  meteorological  instrumentation. 

Three  phase-coherent  channels  at  f ^  =  11.4  GHz,  f2  =  28.8  GHz,  f^  =  96.1  GHz  provided 
amplitude  data  suitable  to  analyze  atmospheric  loss  and  fade  mechanisms  and  evaluate 
■frequency-dependent  responses.  Absolute  calibration  of  the  received  signal  levels 
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C*/f|  s  104  <V>2/<T>1, 


(21) 


and  humidity  variance  dominates  the  refractivity  structure.  The  resulting  c*  values 

are  large  ("moist"  fluctuations).  The  structure  constant  is  highest  in  layers  near* 

•12  -17 

est  the  ground  where  millimeter-wave  values  for  c*  can  range  between  10  and  10 
-2/3  n 

B 

The  refractivity  N  given  by  (3)  is  complex.  Terms  should  ba  added  to  (17)  to 
account  for  contributions  from  the  variance  of  absorption  and  the  covariance  of  re¬ 
fraction  and  absorption  (Ott  and  Thompson,  1978;  Hill  and  Clifford,  1981).  Usually, 
humidity  variations  dominate  the  structure  constant  c*  via  refraction  ,  and  the  con¬ 
tribution  of  absorption  to  fluctuations  will  be  difficult  to  observe  in  the  milli¬ 
meter-wavelength  range,  unless  specific  values  of  about  5  dB/km  are  exceeded  (Filho 
et  al. ,  1983). 

When  a  wind  with  average  speed  u^  (m/s)  transverse  to  the  propagation  path  car¬ 
ries  the  spatial  structure  of  refractivity  fluctuations  through  the  propagation 
path,  a  modulation  spectrum  of  scintillation  frequencies  u»  (Hz)  Is  imposed  upon  the 
propagating  wave.  Between  inner  and  outer  scales  of  turbulence  a  spatial  spectrum 
exists  which  has  a  theoretical  (Kolmogorov)  foundation.  Turbulence  structures 
larger  than  the  size  of  the  1st  Fresnel  zone  produce  signal  scintillations  with  con¬ 
stant  spectral  power  density  W  fdB2/Hz).  These  occur  in  the  lower  frequency  portion 

°  -8/3 

of  the  spectrum.  The  power  density  falls  off  rapidly  (w  )  with  Increasing  fre- 

-8/3 

quencies  u>.  The  spectral  asymptotes  WQ  and  t.  cross  at  the  "corner"  frequency 

u»0  ?  u^/^Xd  =  0.058  u^/T/cT  Hz.  (22) 

The  spectral  power  density  for  the  range  w  »  w0  is  proportional  to  (Cole  et  al., 
1978) 


V(w)  «  f*  c*  d  £j/3  dBVHz.  (23) 

Structure  sizes  larger  than  tQ  are  ill-dafined  in  form.  Turbulent  phenomena 
merge  with  large-scale  variations  in  the  atmosphere.  A  power  spectrum  W(w)  should 
reveal  the  processes  discussed  and  allow  the  deduction  of  a  relative  refractive  struc 
ture  constant  c*  «  WQ(w)  and  w*nd  speed  u  via  (17)  and  (22). 

2.3  Multipath  Fading 

In  some  regions  of  the  world,  frequency- selective  fading  phenomena  from  atmos¬ 
pheric  refractive  layers  are  fairly  common.  They  are  associated  with  stable 
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The  propagation  experiment  described  in  Section  3  yields  time  series  s(t)  of 
log-amplitude  behavior  at  three  frequencies  f.  Under  atmospheric  conditions  where 
(16)  is  valid,  the  variance  of  log- amplitude  scintillations  is  expected  to  be  propor¬ 
tional  to  (Vogel  et  al.,  1984) 

ofi  «  f7/6d11/6c2  dB2  (17) 

a  n 

Alf'O,  it  is  assumed  that  the  transmitting  and  receiving  antenna  sizes  are  small  rela¬ 
tive  to  an  approximate  1st  Fresnel  zone  clearance  Jxd,  and  in  turn,  the  Fresnel  zone 
is  small  compared  with  the  outer  turbulence  scale  The  radio  wavelength  X  fol¬ 
lows  from 

X  =  299. 8/f  mm.  (18) 

Since  both  air  density  and  absolute  humidity  v  affect  the  refractivity  N,  the 
structure  parameter  c*  can  be  expressed  in  terms  of  analogous  parameters  for  temper¬ 
ature,  humidity,  and  barometric  pressure.  Neglecting  pressure  variations, 
c*  is  given  by 


CS  =  Af 


7T&  *  A$ 


*  2ATAv  <T><v> 


where  the  sensitivity  coefficients  Aj  y  are  calculated  from  the  refractivity  NQ  (Kill 
and  Clifford,  1981).  The  last  term  in  (19)  can  be  approximated  by 


_  Tv 
<T><v> 


9tl 


The  positive  sign  usually  applies  during  the  day  and  the  negative  sign  during  the 
night,  due  to  a  change  in  the  direction  of  the  temperature  gradient  near  the  ground. 
Under  some  dry  conditions  (Hebe  and  Hopponen,  1977),  (20)  can  substantially  cancel 
the  contribution  of  the  other  two  terms  and  result  in  a  situation  for  which  the 
value  of  c*  is  small  ("dry"  scintillations).  Frequently,  the  ground  evaporates 
water  vapor.  In  these  instances,  humidity  fluctuations  will  surpass,  in  relative 
terms,  temperature  fluctuations  by  far,  say 
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Figure  2.  Predicted  specific  attenuation  a  and  dispersive  delay  £  for  moist 
air  (0>IOO*RH)  at  sea  level  including  simulated  fog  conditions 
(w  »  0.1)  with  roughly  300  m  visibility.  Five  atmosoheric  mil  1 i > 
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Figure  1.  Predicted  (MPM)  specific  attenuation  a  (a)  and  dispersive  delay  3  (b) 
up  to  1000  GH:  for  moist  sea  level  air  (#1  and  3)  and  pure  water  vapor 
(#2  and  4)  each  for  two  conditions  v  •  1  and  30.3  (100'RH)  g/nr  at  T  »  30*C 
(Note  the  linear-log  scales  used  to  display  delay  6  ) 
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In  s unwary,  (1)  to  (IS)  constitute  the  physical  HPH,  which  Is  capable  of  pre- 
dieting  millimeter  wave  properties  of  »oist  air  for  given  sets  of  P-T-v(RH)-w  distri- 

•4  c 

buttons  along  a  radio  path.  A  full  range  of  HPH  coverage  (i.e. ,  a(f)  =  10  to  10 
dB/km,  p(f)  =  ±10*  ps/km,  f  *  1  to  1000  GHz)  is  exhibited  in  Figure  1  for  two  cases 
[v  =  1  and  30.3(100%RH)  g/mJ ) ,  each  for  pure  water  vapor  (H^0  self-broadening)  and 
moist  sea  level  air.  Examples  in  Figure  2  show  the  humidity  dependence  of  specific 
attenuation  a  and  dispersive  delay  0  for  sea  level  conditions  (P  =  101.3  kPa,  RH  -  0 
to  100S,  w  =  0.1  g/m3,  T  =  5  and  25°C)  and  frequencies  up  to  f  =  350  GHz.  Five  at¬ 
mospheric  transmission  windows  are  Indicated  by  W1  to  W5. 

Starting  from  a  given  choice  of  mean  atmospheric  conditions,  it  is  the  fluctua¬ 
tions  in  P-T-v-w  that  modulate  the  average  moist  air  refractivity  N  t  AN(t,x)  in 
space  ana  time  which,  in  turn,  influences  the  propagating  radio  wave.  To  evaluate 
the  resulting  propagation  effects,  it  is  necessary  to  investigate  the  stochastic 
functions  that  govern  atmospheric  fluctuations.  The  practical  effects  of  refractiv¬ 
ity  N  and  its  variance  a  are  to  produce  phenomena  that  may  be  grouped  into  those 
resulting  from  small-scale  turbulence  structures  (see  Section  2.2)  and  those  from 
large-scale  structures  in  the  atirasphere,  such  as  diurnal  and  weather  changes,  over¬ 
all  temperature  and  water  vapor  distributions  including  inversions,  etc.  The  latter 
manifest  themselves  in  effects  such  as  ray  bending,  ducting,  propagation  delay  dis¬ 
persion,  and  multipath  interference  (see  Section  2.3). 

2.2  Atmospheric  Turbulence 

The  theory  of  wave  propagation  in  a  turbulent  medium,  as  described  by  Tatarski 
(1961),  was  extended  to  millimeter  waves  by  Clifford  and  Strohbehn  (1970),  Ishimaru 
(1972),  Ott  and  Thompson  (1978),  and  Hill  and  Clifford  (1981),  among  others.  Sub¬ 
millimeter  wave  turbulence  effects  have  been  treated  recently  in  a  review  by  Bohlander 
et  al.  (1985).  Under  ideal  assumptions  that  turbulent  eddies  are  distributed  within 
a  propagation  path  in  an  isotropic  and  homogeneous  manner,  the  structure  function  at 
separation  Ax  can  be  written  as 

o  =  <aN2>  =  c2  ax2/3  1012,  (16) 

where  <  >  denotes  ensemble  averages,  usually  replaced  by  time  averages.  The  separa¬ 
tion  Ax  of  two  measurement  points  has  to  be  between  an  inner  scale  *.  (<  0.01  m)  and 

2  -2/3  % 

an  outer  scale  *o(>  5  ra),  and  cR  in  units  of  m  is  called  the  refractive  struc¬ 
ture  constant.  This  parameter  is  the  key  in  relating  turbulence  to  propagation  ef¬ 
fects  such  as  amplitude  scintillations,  phase  front  distortions,  and  variations  in 
angle-of-arrival. 


can  bt  used  Instead  of  (14a),  based  on  data  reported  by  Simpson  et  al.  (1979). 

The  dielectric  data  c',e"  of  bulk  water  are  calculated  with  the  Debye  model, 
which  is  valid  for  f  <  300  GHz  as  reported  by  Chang  and  Wilheit  (1979): 

c"  *  (185  -  113/fli)ft|/[l  ♦  (ftw)2)  (15a) 

and 

c*  «  4.9  ♦  (185  -  113/8)/[l  ♦  (ft^)2],  (15b) 

where  t  =  4.17  *  10~^  8  exp(7.13  8)  ns.  Table  3  lists  numerical  examples  of 
aw  *  0. 182f  N^  (dB/ka)  for  two  temperatures  (0°C  and  25°C)  at  selected  frequencies. 

A  strong  temperature  j  pendence  can  be  noticed. 

Haze  conditions  are  related  to  vapor- to*droplet  conversion  processes  reversible 
(swelling/shrinking)  with  relative  humidity  for  values  below  RH  =  ICO  percent.  The 
simple  growth  function  (5)  provides  an  approximation  of  the  RH-reversible  conversion 
up  to  RH  a  96  percent.  Typical  dry  air  mass  loadings  of  hygroscopic  aerosol  in 
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ground  level  a*r  are  below  wQ  3  10  g/m  .  A  practical  growth  function  covering  the 
range  above  RH  *  %  percent,  including  the  irreversible  condition  of  supersatura¬ 
tion,  RH>100  percent,  is  lacking.  Such  a  function  is  needed  to  explain  the  growth 
of  suspended  droplet  concentrations  to  values  w  =  10"3  to  1  g/m3  as  observed  under 
various  fog  conditions.  At  present,  reported  values  for  w  are  used  in  (14)  and  the 
contributions  are  added  to  respective  calculations  (7)  for  saturated  (100XRH)  air. 

Table  3. 

Haze  and  Fog  Attenuation  ow  and  Delay  Pw  Normalized  to  w  *  1  g/m3 
TEMPERATURE  FREQUENCY  f  (GHz) 


T 

1 

10 

30 

100 

200 

300 

400 

600 

800 

1000 

1 

#C 

0 

.0010 

.097 

.82 

5.4 

(dB/ksi) 

9.3  10.8 

13 

18 

23 

29 

25 

.0005 

.051 

.45 

4.2 

10.8 

15.3 

21 

31 

40 

48 

4  1 

0 

.69 

.32 

.09 

.04 

Pw  fps/k ») 

.04  .04 

0 

0 

0 

0 

•*, 

25 

.62 

.49 

.20 

.06 

.04 

.04 

0 

0 

0 

0 

J 

(14  a.b)  (14c) 


,«ake  a  saall  contribution  at  ground  level  pressures  due  to  the  nonresonant  0^  spec* 
tr»«  below  10  GHz  and  a  pressure* induced  spec true,  effective  above  100  GHz 

(Stankevich,  1974).  A  width  para^ter  for  the  Oebye  spectrua  of  0^  is  foraulated  in 
accordance  with  (10)  to  be  s  5-6  *  10_3(p  ♦  ).le)8°*8  (GHz)  (Rosenkranz,  1982). 

0  *4 

The  continuue  coefficients  are  a0  =  3.37  x  10  ,  (Rosenkranz,  1975),  and 

ap  »  1.17  x  10’10  (Oagg  et  al..  1982). 

The  water  vapor  continuua  is  an  eepirical  contribution  derived  froa  fitting  ex- 
periaental  data  in  the  case  of  N"  ana  theoretical  data  in  the  case  of  N^,  which  led 
to 

«;(f)  *  [bfp  ♦  bee©3]f  e  825  (13a) 

and 

N;(f)  a  V  02‘4f2‘05,  (13b) 

where  bQ  *  6.47  x  10‘8  (R.  Hill,  personal  conauni cat ions,  1984), 
bf  »  1.40  x  10-6  and  b#  ■  5.41  x  lo'5  (Liebe,  1984). 

2.1.3  Hyorosol  Continuua  (Haze  and  Fog) 

Suspended  water  droplets  (hydrosols)  in  haze  and  fog  (or  clouds)  are  ailliaeter 
wave  absorbers.  Their  size  range  of  radii  is  below  50  pa,  which  allows  the  Rayleigh 
approxiaalion  to  Hie  scattering  theory  to  be  used  for  calculating  refractivity  con¬ 
tributions  Nw  to  (7)  in  the  fora  (liebe,  1981;  Falcone  et  al.,  1982) 

N“(f)  ’  4.50w/c"(I  ♦  C2)  (14a) 

and 

fT(f)  ^  2.4  *  10’3w  c’,  (14b) 

where  C  *  (2  ♦  c’)/c";  and  c\c"  are  real  and  imaginary  parts  of  the  dielectric  con¬ 
stant  for  water,  for  frequencies  above  300  GHz,  the  following  approximation 


ifl 


The  line  parameters  are  calculated  according  to  the  scheme  below: 

Syabol 

02  lines  in  air 

H^O  lines  in  air 

S  .  kHz 

«jP03  10’6  exp[a2(l  -  0)) 

bjee3,5  exp£b2(l  -  0)) 

(9) 

Y  .  6Hz 

(0.8-a.)  - 

a3(p0  4  ♦  i.le6)10  J 

b3(pe0’8  ♦  4.80e6)10‘3 

(10) 

6 

®£  -3 

a5p0  *  10  J 

0 

(11) 

Line  center  frequencies  wQ  and  the  spectroscopic  coefficients  to  a^,  bj  to  b3  for 
strength  S,  width  y,  and  overlap  correction  £  are  listed  in  Table  2. 

Standard  line  shapes  FM(f),  including  the  Van  Vleck-Weisskopf  function  (8), 
predict  in  frequency  regions  of  local  line  dominance  about  the  sane  results  for 
NM(f)  as  long  as  f“(f)  values  are  larger  than  0.1  percent  the  peak  value  at  f  *  vQ. 
Wing  contributions  that  are  less  than  0.1  percent  depend  upon  the  chosen  line  shape 
function.  Far-wing  contributions  (<  0.001  percent)  differ  by  as  much  as  1:80 
(Hebe,  1S84).  No  line  shape  has  been  confiraed  that  predicts  absorption  intensi¬ 


ties  over  ranges  (10*3  to  <10*6)  of  F"(v  ),  as  required  for  the  water  vapor  spec- 


trua.  Expected  far-wing  contributions  froa  st;ong  (>  10  d3/ka)  infrared  iines  are 
accounted  for  summarily  by  eapirica!  correction. 


2.1.2  Continuua  Soectrua  for  Air 

Continuua  spectra  in  (7)  identify  dry  air  and  water  vapor  teras  N  ♦  N  and 

P  * 

aust  be  added  to  the  selected  group  of  local  0^  and  resonance  lines  (Table  2) 
described  by  (8)  in  order  to  correctly  predict  ataospherlc  ailliaeter  wave  attenua¬ 
tion  In  window  ranges  between  lines.  Continuua  absorption  increases  monotonically 
with  frequency. 

The  dry  air  continuua 


Np(f )  =  (?«0(T0[1  *  (f/Y0>?JU  ♦  (f/M)2])'1  ♦  P  ap  jf  p  82 


(lie) 


Wp( f )  *  e0{(l  ♦  (f/y0)?f 1  *  UP  0 1 


(12b) 


3 


if 


i,i:v:).U^AOiWJt^L'WA!?l?Wi>»?il<!i^»li.i|l.'i!1 


N"(f)  =  £  (SF“),  ♦  N"  ♦  £  <SF"),  ♦  N"  ♦  N" 
1=1  ’Pi  tew 


H'(f)  =  £  (SF‘) 

fa 


<*vi! 


(SF*)t  ♦  n;  ♦  n; 


dry  air  ♦  water  vapor  ♦  susp.  droplets. 


In  the  line-by-line  calculations,  S  Is  a  line  strength  in  kHz  and  F*  and  F“  are  real 
and  laaginary  parts  of  a  shape  function  in  GHz"*.  Pressure  broadening  of  the  domi¬ 


nant  lines  leads  to  two  different  types  of  frequency  responses,  namely  sharp  reson¬ 
ance  lines  and  continuum  spectra.  Both  are  treated  in  the  following  for  ataospheric 
conditions  up  to  heights  of  h  *  30  ka.  Self-  and  foreign-gas*broadening  influences 
have  to  be  taken  into  account.  Trace  gas  spectra  are  neglected  since  their  abundance 
is  too  saall  to  markedly  affect  tropospheric  propagation. 


2.1.1  Local  Line  Absorption  and  Dispersion 

The  Van  Vleck-Weisskopf  function  was  modified  by  Rosenkranz  (1975)  to  describe, 
to  first  order,  line  overlap  effects.  This  leads  to  local  absorption  and  dispersion 
line  profiles 
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Figure  3.  Terrain  profile  of  27.2-km  Boulder  path  with  two  ray  paths  for  normal 
(k=l . 34)  and  subrefracti ve  (k=0. 29)  propagation  (Espeland  et  a  1 . ,  1984; 
Stephansen,  1981): 

RX  =  receiver  site  in  center  of  Boulder 
TX  =  transmitter  site  north  of  town 


Full  details  of  the  transmitter  and  receiver  configurations  have  been  given  by 
Espeland  et  al.  (1984).  All  salient  parameters  that  are  needed  to  make  a  perfor¬ 
mance  assessment  (24)-(29)  for  each  of  the  three  channels  are  summarized  in 
Table  4.  The  cw  carriers  are  phase-locked  to  a  100  MHz  crystal  oscillator.  Trans¬ 
mitters  and  receivers  both  are  located  inside  protective,  temperature-controlled  (39° 
±  1  °C)  structures.-  Link  gain  stability  is  maintained  at  ±  0.5  dB  over  an  ambient 
temperature  range  from  -25°C  to  +  35°C. 

Meteorological  data  are  recorded  only  at  the  receiving  terminal  which  is,  admit¬ 
tedly,  insufficient  to  capture  atmospheric  activity  over  the  entire  path.  The  mete¬ 
orological  instrumentation  measures  in  situ  (RX)  microwave  (9.375  GHz)  refractivity 
Nq  (6), ‘ambient  air  temperature  T,  and  barometric  pressure  P.  From  these  parameters, 
water  vapor  pressure  e  and  relative  humidity  RH  are  calculated  using  equation  (6), 
(30).  The  following  uncertainties  for  the  sensors  have  been  estimated: 

ANq  =  ±  0.5  ppm,  AT  =  ±  0.5°C,  and  AP  =  ±  0.1  (Ae  =  ±  0.03)  kPa. 

3.2  Data  Acquisition 

All  data  logging  and  recording  operations  are  controlled  by  a  mini-computer  as 
shown  schematically  in  Figure  4.  Input  signals  are  sampled  by  a  scanning  A/D  conver¬ 
ter  (5  mV  accuracy;  12  bit  resolution;  55,000  readings  per  second).  The  amplitudes 
S1  2  3  of  the  cw  carriers  can  be  sampled  at  rates  up  to  10  per  second.  To  reduce 
the  amount  of  data  to  be  stored,  daily  records  were  produced  at  rates  of  1  per 
second,  including  the  meteorological  readings  of  Nq,  P,  T  (and  rain  rate  R,  mm/h). 
Data  initially  is  stored  in  hourly  segments  on  hard  discs  and  periodically  (about 
every  two  weeks)  transferred  to  magnetic  tape  for  archiving.  Time  series  of  eleven 
channels  (s^,  s^,  s^,  wide-band  response,  log  (BER),  e,  RH,  Nq,  T,  P,  R)  for  the 
current  day  are  displayed  on  the  computer  terminal  and  copiecTonce  a  day.  Examples 
from  the  underlined  channels  have  been  selected  for  analysis. 

3.3  Calibration  Procedure  (0.1-km  Path) 

An  overall  system  stability  for  amplitude  detection  of  better  than  ±0.5  dB  was 
maintained  for  all  of  the  three  channels  of  the  Boulder  LOS  experiment  during  the 
measurement  period  (August  1983  -  June  1984).  Calibrated  mean  signal  levels  per¬ 
mitted  studies  of  water  vapor  attenuation  at  28.8  and  96.1  GHz.  Clear  air  amplitude 
scintillations  have  been  recorded  continuously  with  a  sampling  rate  of  1  Hz.  On  a 
few  occasions  the  rate  was  set  to  10  Hz  (maximum). 
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Table  4. 

Performance  Parameters  for  Boulder  LOS  Link 
(Espeland  et  al.,  1984) 

CHANNEL 
2 


)  A,  dB 


E’  K 

,f,  MHz 

...  dBm 


Equation 


142.3 

x  = 

11.4 

27.2 

75 

18.75 

0.46 

32.0 

(4  deg) 

0.60 

0.46 

32.0 

(4  deg) 

-2.0 

-3.5 

-8.5 

68.75 

0.33 
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-74.5 

X 
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0.0015 
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60 

X 

27.7  x 
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-1.0 
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10.3  Figure  6 

10.1  (0.371  dB/km)  Figure  13 
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During  the  experiments  the  96.1  GHz  system  was  carefully  calibrated  four  times. 
This  was  accomplished  by  shortening  the  propagation  path  length  to  d(>  d  1  =  0.105 
km:  free  space  loss  Lq  is  reduced  by  48.3  dB  (from  160.8  to  112.5  dB)  and  atmos¬ 
pheric  loss  A  <  0.05  dB  =  0  is  eliminated.  The  LOS  system  gain  was  compensated  at 
the  receiver  by  inserting  an  accurately  (NBS)  calibrated  attenuator  (58.3  ±0.1  dB) 
to  restore  the  originally  received  level  to  s3  =  -10.0  dB.  In  the  short  path  con¬ 
figuration,  experiments  were  performed  to  evaluate  the  96.1  GHz  receiver  perform¬ 
ance.  The  receiver  noise  power  at  the  rf  mixer  was  determined  to  be  P^  =  -64.0  dBm 
(see  link  budget  in  Table  4). 

When  the  original  path  distance  d  =  27.2  km  was  restored,  only  atmospheric  loss 
A  and  pointing  error  loss  were  left  undetermined.  Repetitive  scans  in  elevation  and 
azimuth  at  both  the  transmitting  and  receiving  terminals  assured  main-beam  align¬ 
ment.  Signal  loss  from  pointing  errors  was  reduced  to  less  than  -0.5  dB.  This 
level  of  accuracy  is  accomplished  by  selecting  optimum  times  for  pointing  (during 
periods  of  minimum  scintillation)  and  by  signal  averaging.  The  overall  uncertainty 
of  the  link  calibration  was  estimated  to  be  smaller  than  ±0.03  dB/km.  For  example, 
an  experiment  performed  at  1315  MST  on  9  August  1983  gave  a  differential  path  loss 
of  58.3  dB,  by  coincidence  the  same  value  as  the  attenuator  used  for  the  calibration 
path.  The  free  space  loss  difference  is  48.3  dB,  hence  the  atmospheric  loss 
A  =  10.0  dB.  The  meteorological  conditions  at  the  receiver  site  were  P  =  83.4  kPa, 

T  =  27°C  and  Nq  =  260  ppm.  Equation  (6)  can  be  rearranged  to  read 

v  =  (Nq  -  2.588  P  0)/(5. 768  0  -  0.028)  g/m3 ,  (30) 

which  yields  v  =  7.69  g/m3  and  leads  to  excellent  agreement  with  the  predicted  value 
A  =  10.3  dB  (see  Figure  6  and  Table  4).  Detection  of  path  loss  A  at  28.8  GHz  was 
possible  by  a  relative  calibration  to  s^  =  0  dB  from  a  "best"  fit  of  the  experimental 
data  to  predicted  values  (see  Figure  13). 

4.  OBSERVATIONS,  PREDICTIONS,  AND  RESULTS 

In  clear-air  propagation  of  millimeter  waves,  weather-related  changes  of  atmosphe¬ 
ric  refractivity  (N)  are  mainly  responsible  for  level  changes  and  short-term  fluctua¬ 
tions  of  received  signal  amplitudes  (s^).  The  propagation  effects  of  attenuation, 
scintillation,  and  delay  dispersion  within  the  signal  bandwidth  reduce  the  fade  mar¬ 
gin  (F)  of  a  system.  The  program  MPM,  described  in  Section  2,  is  a  valuable  aid  to 
support  planning  decisions  for  anticipated  systems  and  to  make  performance  assess¬ 
ments.  On  a  more  basic  scale,  measurements  of  a  and  p  reveal  information  on  molecular 


interactions,  and  RSL  scintillations  offer  insights  into  atmospheric  dynamics. . 

Experiments  have  been  conducted,  including  the  one  described  in  Section  3,  pro¬ 
ducing  data  that  are  useful  to  verify  various  aspects  of  the  highly  structured  mil¬ 
limeter-wave  refractivity  N  (i.e. ,  program  MPM).  In  this  section,  results  are  pre¬ 
sented  regarding  frequency-humidity-temperature  dependences  of  attenuation  and,  for 
one  case,  of  associated  delay  dispersion.  From  the  Boulder  LOS  experiment,  the 
daily  records  depicted  in  Figure  5  are  presentations  of  minimum/maximum  and  strong¬ 
est  anomalous  scintillation  behaviors,  selected  from  almost  continuous  records 
covering  an  11  month  period.  The  spectral  density  W(u>)  of  amplitude  scintillations 
at  11.4,  28.8,  and  96.1  GHz  during  these  events  is  calculated. 

4.1  Received  Amplitude  Levels 

Time  series  s^  2  j(t)  have  been  recorded  and  particular  events  are  analyzed 
for  the  behavior  of  mean  s  and  variance 

og  =  (n  -  1)  1  (Asn/s)2  dB2, 

with  n  the  number  of  samples.  Daily  summary  records  containing  selected  event  seg¬ 
ments  are  depicted  in  Figure  5  for  four  cases:  (a)  a  relatively  quiet  period  ("dry" 
scintillations),  (b)  a  relatively  active  period  ("moist"  scintillations),  and  (c)  a 
stretch  of  rapid  and  severe  up-and-down  fades  (possible  multipath  effects).  In  ad¬ 
dition,  5d  displays  a  rain  event  to  demonstrate  its  severe  influence--but  precipita¬ 
tion  will  not  be  a  topic  in  any  discussion  to  follow.  A  statistical  analysis  usu¬ 
ally  reveals  a  log-normal  cumulative  probability  distribution  for  periods  (~  10  min) 
of  stationary  conditions  (e.g..  Cole  et  al.,  1978). 

The  power  spectra  of  the  amplitude  scintillations  were  computed  using  the  met¬ 
hod  of  averaged  modified  periodograms  (Oppenheim  and  Schafer,  1975).  The  Blackman 
window  (Harris,  1978)  was  used  and  normalized  for  a  coherent  gain  of  unity  (x  2.38). 
Modified  periodograms  were  computed  for  segments  of  256  points  (sampled  at  a  1  Hz 
rate)  to  estimate  the  high  frequency  portion  of  the  spectrum.  The  data  for  the  time 
period  of  interest  were  then  averaged  in  blocks  of  256  points  and  segments  of  16 
points  used  to  estimate  the  low  frequency  portion  of  the  spectrum. 

For  purposes  of  clarity  in  the  graphical  presentation  of  the  ratios  of  the  spec 
tra,  the  averaged  periodograms  were  smoothed  so  that  the  points  would  be  equally 
spaced  in  frequency  on  the  logarithmic  scale  before  the  ratios  were  calculated. 

This  was  done  by  block  averaging,  beginning  at  the  spectral  point  lowest  in  fre¬ 
quency  and  increasing  the  number  of  points  averaged  as  the  frequency  increased. 
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Figure  5.  continued  (c) 
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Received  CW  Signal  Levels  (path-average) 

(d)  Rain  event  (1800-1900,  06/06/84): 


Specifically,  the  number  of  points  in  each  block  was  doubled  beginning  with  a  block 
of  one  point  for  the  lowest  frequency. 

Amplitude  scintillations  are  affected  by  temperature  and  humidity  fluctuations 
along  and  wind  across  the  radio  path.  To  check  how  well  the  assumptions  of  turbu¬ 
lence  theory  [see  above  (16)]  were  met,  data  could  be  drawn  from  an  extensive  net¬ 
work  of  meteorological  observation  stations.  In  situ  climatology  within  5  km  proxi¬ 
mity  of  the  test  path  was  recorded  by  five  different  groups: 


Acronym 

Organization 

#  of  stations 
available 

LOS 

Boulder  test  link,  receiver  site 

NTIA-ITS.S 

1  (see  Fig.  5) 

BAO 

Boulder  atmospheric  observatory 

N0AA-WPL 

1  (10  m  level) 

BMW 

Boulder  wind  observatory 

N0AA-WPL 

4 

PROFS 

surface  mesonetwork 

N0AA-ERL 

3 

MAYPOLE 

surface  mesonetwork 

NCAR-MRS 

6 

and  made  available  to  ITS.  Special  thanks  for  exemplary  cooperation  is  in  order 
for  Dr.'s  Chandran  Kaimal  (BAO),  Alfred  Bedard  (BMW),  Peter  Mandics  (PROFS),  and 
Ms.  Cindy  Mueller  (MAYPOLE). 


4.2  Predicted  Responses 

Clear  air  transmission  effects  formulated  by  (l)-(3)  represent  ultimate  perfor¬ 
mance  limitations  to  millimeter-wave  applications.  The  MPM  program  was  developed  to 
be  valid  for  frequencies  up  to  1000  GHz.  Computer-efficient  operation  is  assured  by 
using  only  a  limited  local  line  data  base  (Table  2)  of  the  major  absorbers  H^O  and 
O2,  supplemented  by  continuum  spectra  (12)-(14)  including  a  suspended  droplet  (haze, 
fog)  term.  Humidity  response  of  specific  attenuation  a(v)  was  predicted  for  the 
Boulder  LOS  link  as  displayed  in  Figure  6.  Useful  data  to  confirm  predictions  can 
be  expected  only  at  96.1  GHz  since  the  gain  stability  of  the  test  link  is  on  the 
order  of  ±0.05  dB/km. 

The  combined  influence  of  both  absolute  (v)  and  relative  (RH)  humidity  upon 
attenuation  is  made  clear  by  analyzing  the  temperature  dependence  of  a  water  .apor 
attenuation  slope 

g  =  [cr( v)  -  a(v=0)]/v  (dB/km)/(g/m3 ) . 

Results  of  g(T)  are  presented  in  Figure  7  for  three  cases:  (a)  for  frequencies  of 
the  Boulder  LOS  link,  (b)  for  two  water  vapor  line  peaks,  and  (c)  for  centers  of 
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Figure  7.  continued... 

b)  H^O  line  center  frequencies  (2x) 

c)  atmospheric  transmission  window  frequencies  (4x) 
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four  atmospheric  windows  Wl  to  W4  (see  Figure  2).  Corresponding  numerical  values  of 
specific  attenuation  a(RH)  in  Table  5  span  a  wide  range  at  sea  level  conditions. 

MPM  predictions  have  to  cover  reliably  values  between  0.01  and  >  140  dB/km. 

Attenuation  a due  to  suspended  droplets  (Figure  2,  curve  #9)  is  added  to  satur¬ 
ated  air  conditions  that  attenuate  at  rates  cts(100%RH).  The  selected  droplet  concen¬ 
tration  of  w  =  0.1  g/mJ  is  representative  of  a  fog  with  about  0.3  km  optical  visibil¬ 
ity.  A  hydrosol  factor 

xw  3  (V°s)(vs/w) 

is  introduced  in  Table  5  to  demonstrate,  particularly  in  the  window  ranges,  a  strong 
increase  in  attenuation  per  unit  of  H^O  absorbers  when  vapor  molecules  convert  into 
droplets.  Typically,  w  S  0.01  for  haze  and  w  >  0.1  to  1  for  fog  conditions. 

A  strong,  nonlinear  influence  of  atmospheric  humidity  RH(T),  v(RM),  and  w(RH) 
on  millimeter-wave  attenuation  is  predicted  and  elucidated  by  selected  examples  (see 
Figures  6,  7  and  Table  S).  To  what  extent  experimental  evidence  is  supported  by  MPM 
predictions  will  be  investigated  in  the  remainder  of  this  report. 

4.3  Model-Versus-Reported  Data 

Corroborative  experimental  data  of  sufficient  quality  to  scrutinize  predictions 
are  scarce.  Reliability,  precision,  and  limited  scope  of  supporting  meteorological 
data  often  compromise  the  accuracy  of  results  deduced  from  field  observations.  Gen¬ 
erally,  laboratory  experiments  provide  more  precise  tests  by  simulating  controlled 
electromagnetic  and  atmospheric  conditions  crucial  to  model  validations.  In  this 
manner,  contributions  from  the  22  and  183  GHz  water  vapor  lines  and  48  to  70,  119 
GHz  oxygen  lines  have  been  evaluated  as  summarized  by  liebe  (1981,  1983,  1985). 
Confirmation  Is  lacking  in  window  ranges  W2  to  W5  (see  Figure  2)  that  are  sensitive 
to  the  water  vapor  continuum  (13).  To  test  the  MPH  program,  reported  data  sets  on 
humidity-dependent  propagation  effects  were  chosen.  The  data  have  been  gathered 
recently  at  different  locations  under  various  conditions.  An  overview  of  the  experi¬ 
ments  is  given  In  Table  6.  Reported  data  are  generally  available  in  graphic  form 
and  had  to  be  digitized  by  us  for  a  comparison  with  model  predictions. 

The  first  example,  Figure  8,  is  unique--not  only  are  attanuation  rates  e(v) 
reported  at  the  frequencies  fj  =  81.84  and  ■  245.52  GHz,  but  also  included  are 
data  on  the  associated  differential  dispersion  (7b),  (2b) 

A*(v)  *  ap/3.336  «  Nf(f2)  -  N’(fj). 
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Table  5. 

Predicted  Specific  Attenuation  a  (dB/km)  for  a  Total  Pressure  P  *  101.3  kPa  and 
Hydrosol  Attenuation  Factor  x  at  Various  Relative  Hunidities  RH  (X)  and  Teaperatures 
T(K)  for  Selected  Milliaeter-Bave  Frequencies 


Parameters 

f  T  vs 

Hydrosol 

Factor 

*w 

10CX  RH 

Specific  Attenuation  or  of  Air 

75%  RH  50X  RH  25X  RH  OX  RH 

GHz 

K 

g/a3 

<ffl/  ka 

22.2 

310 

43.46 

8.4 

1.03 

0.78 

0.52 

0.27 

BSSl 

25.49 

0.62 

0.46 

0.16 

(vj 

290 

14.31 

13 

0.35 

0.27 

0.18 

0.10 

0.013 

0 

280 

7.65 

17 

0.19 

0.15 

0.10 

0.06 

0.014 

270 

3.87 

22 

0.11 

0.08 

0.06 

0.04 

0.016 

260 

1.85 

29 

0.06 

0.05 

0.04 

0.03 

0.017 

35.0 

310 

26 

0.76 

0.50 

0.29 

0.13 

0.026 

300 

41 

0.38 

0.27 

0.17 

0.09 

0.028 

(Wl) 

290 

m 

60 

0.20 

0.15 

0.10 

0.06 

0.031 

280 

84 

0.12 

0.09 

0.07 

0.05 

0.034 

270 

114 

0.08 

0.07 

0.06 

0.05 

0.038 

260 

135 

0.06 

0.06 

0.05 

0.05 

0.042 

95.0 

310 

28 

4.56 

2.89 

1.58 

0.63 

0.036 

300 

44 

2.18 

1.44 

0.83 

0.37 

0.040 

(W2) 

290 

m 

62 

1.05 

0.73 

0.45 

0.22 

0.044 

280 

79 

0.53 

0.38 

0.26 

0.14 

0.048 

270 

90 

0.28 

0.21 

0.16 

0.10 

0.053 

260 

89 

0.16 

0.14 

0.11 

0.08 

0.058 

140.0 

310 

24 

10.21 

6.48 

3.54 

1.39 

0.019 

300 

35 

4.88 

3.21 

1.84 

0.78 

0.021 

<W3) 

290 

« 

46 

2.34 

1.60 

0.97 

0.44 

0.023 

280 

53 

1.13 

0.80 

0.51 

0.25 

0.025 

270 

54 

0.56 

0.41 

0.27 

0.15 

0.027 

260 

48 

0.28 

0.21 

0.15 

0.09 

0.029 

183.3 

310 

3.0 

143.08 

109.98 

75.43 

38.96 

0.014 

300 

2.9 

91.32 

69.60 

47.22 

24.07 

0.016 

(O 

290 

• 

2.7 

54.94 

41.59 

28.01 

14.16 

0.017 

0 

280 

2.5 

31.10 

23.44 

15.71 

7.91 

0.018 

270 

2.1 

1G.54 

12.44 

8.32 

4.18 

0.019 

260 

1.7 

8.24 

6.19 

4.14 

2.03 

0.020 

220.0  310 

18 

26.23 

16,73 

9.19 

3.62 

0.016 

300 

23 

12.64 

8.35 

4.01 

2.04 

0.018 

<W)  290 

27 

6.10 

4.18 

2.52 

1.14 

0.019 

230 

28 

2.97 

2.10 

1.32 

0.63 

0.021 

270 

25 

1.45 

1.06 

0.69 

0.34 

0.022 

260 

20 

0.70 

0.53 

0.35 

0. 19 

0.023 

Curve  10  (Fig.  2, 
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Table  6. 


Overview  of  Experimental  Results  From  Horizontal  Line-of-Sight  (LOS)  Paths,  a 
Mountain  Peak  Zenith  (ZEN)  Path,  and  Laboratory  (LAB)  Experiments,  All  Used  As 
Corroborative  Oata  for  the  Propagation  Program  MPM 


Test  Frequency  f  Figure  Path  Type 


GHz 

81.84 

8 

LOS 

245.52 

8 

LOS 

337 

9 

LOS 

94 

10 

LOS 

192-260 

11 

LOS 

335-416 

12a 

LOS 

330-430 

12b 

LAB 

138.2 

(Eq.  13a) 

LAB 

96.1 

13.  6 

LOS 

2.5-10.0 

(Table  7) 

ZEN 

a  *  Manabe  et  al.t  1984 
b  *  Gasfewski,  1983 
c  *  Buijs  and  Janssen,  1981 
d  •  Fedoseev  and  Koukin,  1984 
e  ■  Furashov  et  «!.,  1984 
f  •  Liebe.  1984 

g  *  L.  Oanese,  personal  communication 


Path  Length  d 

Expon.  y  Reference 

km 

(T*1  law) 
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3. 6(1. 3) 
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Figure  8.  Specific  attenuation  a  and  refractive  dispersion  5N  =  N'(f2)  -  N‘(f,) 
for  humid  air  and  temperatures  every  10°C  ranging  from  G*to  40°C  at  two 
frequencies  fj,?  of  a  Tokyo  propagation  path: 

•  Representative  of  measured  data  clusters  (Manabe  et  al.t  1984) 
—  MPM.  •  RH  =  100* 
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The  curves  <n  Figure  8  were  produced  by  simulating  reported  meteorological  condi¬ 
tions  and  then  superimposing  the  experimental  data.  Clusters  of  original  points 
were  condensed  into  dots  in  the  course  of  the  digitizing  process.  The  gain  stabil¬ 
ity  of  this  LOS  experiment  was  estimated  to  be  better  than  10.3  dB/km  over  a  temper¬ 
ature  range  of  0°  to  40°C.  Good  agreement  is  displayed  between  model  and  experi¬ 
mental  data  on  both  attenuation  and  dispersion.  When  contrasted  with  the  attenua¬ 
tion  rates  o(v),  an  almost  linear  dependence  with  absolute  humidity  v  can  be  ob¬ 
served  for  dispersion  AN(v).  Theoretically,  outside  a  line  center  frequency  range 
vQ  l  By,  refractive  dispersion  N'(f)  is  independent  of  a  particular  shape  function 
F'(f)  which  makes  the  correction  (13b)  for  H^O  line  contributions  above  1  THz  quite 
reliable.  Equation  (13b)  was  formulated  by  fitting  line-by-line  calculations 
(N* (H^O  spectrum)  -  N‘( lines  <  1  THz)]  that  were  made  available  by  R.  Hill  (personal 
communication,  1984).  The  next  example,  Figure  9,  pertains  to  an  LOS  experiment 
performed  at  337  GHz  in  Sandusky,  Ohio,  at  temperatures  from  12* to  25°C.  Figure  9 
shows  one  data  set  o(v),  out  of  18,  selected  over  a  narrow  temperature  range  16°  1 
1#C. 

A  different  type  of  effect  is  often  neglected  in  millimeter  wave  propagation, 
namely  that  depending  on  rela.ive  humidity.  Atmospheric  air  always  contains  a  cer¬ 
tain  mass  concentration  w  of  hygroscopic  aerosol;  for  example,  some  typical  values 


wQ  (RH  <  10%)  »  10  *  (cl##n  a{r)  to  5  x  10  4  (pollute(,  air)  Q/«5. 
w  (RH  3  50%)  »  10  *  to  1  x  10  *  g/m». 

Water  vapor-to-droplet  conversion  follows  roughly  (5)  up  to  96%  RH.  Growth  func¬ 
tions  valid  for  the  range  RH  >  96  to  101%  are  not  easily  available  (e.g. ,  see  Haenel, 
1976).  Oata  of  a(v,RH),  acquired  close  to  the  seashore  in  The  Netherlands  at  both 
94  GHz  and  optical  frequencies  (visible  light)  are  displayed  in  Figure  19.  All  val¬ 
ues  exceeding  a(v)-prcdictions  are  bunched  close  to  RH  a  100%  where  haze  transforms 
into  fog.  Hence,  attenuation  excesses  are  explained  via  (14a)  by  hydrosol  concentra¬ 
tions  up  to  w  3  0.2  g/e1.  Optical  attenuation  (Hie  scattering  theory)  places  the 
boundary  (visibility)  between  ha:e  and  fog  at  about  5  dB/km. 

Water  vapor  attenuation  data  for  frequencies  between  195  and  340  GHz  are  com¬ 
pared  with  model  predictions  in  Figures  11  and  12.  Results  for  the  range  f  *  330  to 
430  GHz  include  field  (8.5°C,  8.4  g/m>,  99%RH)  and  laboratory  (25.5#C,  19  g/m*. 

80%RH)  data  taken  in  the  U.S.S.R.  (Figure  12).  Considering  the  multifarious  diffi¬ 
culties  that  plague  absolute  cal ibratlons,  the  comparisons  of  mode I -vs. -experimental 
data  are  encouraging.  No  anomalous  absorption  features  have  been  uncovered. 
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Figure  9.  Specific  attenuation  a  for  huaid  air  at  T  =  16  ±  2° C  and  a  frequency 
Of  337  GHz: 


♦  Measured  data  (Gasievski,  1934) 
—  MPH.  •  RH  *  100X 
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Figure  12.  Water  vapor  attenuation  rates  o(v)  across  the  atmospheric  window  ranges 
W5  and  W6  at  two  temperatures,  8.5*and  25.5°C: 

♦  Measured  data  (Furashov  et  al.,  1984) 

-  MPM 


■  Vj 


Ww 


Analysis  of  raw  data  a(v,T)  available  from  the  experiments  allows  an  estimate 
of  the  negative  temperature  dependence  0^.  The  exponent  y  was  found  to  vary  between 
3.5  and  6  (Table  5),  which  is  consistent  with  semiempirical  theoretical  results 
(Liebe,  1985). 

Path  attenuation  A  »  Ap  ♦  Ay  and  atmospheric  noise  expressed  by  a  brightness 
temperature  Tg  =  Tp  ♦  Ty  are  closely  related  (e.g. ,  Waters,  1976).  Brightness  temper* 
aturcs  Tg  have  been  measured  at  six  frequencies  between  2.5  and  90  GHz  in  clear 
weather  with  radiometers  looking  toward  zenith  from  a  mountain  peak.  At  hQ  -  3.80  km, 
a  typical  tropospheric  water  vapor  content  V  (mm)  is  reduced  to  about  one-fifth.  In 
fact,  an  average  water  vapor  content  was  determined  to  be 


OB 

■£ 


v(h)dh 


3  ±  1 


All  experiments  were  conducted  with  utmost  care  for  the  purpose  of  measuring  cosmic 
background  radiation  (2.8  K). 

The  MPM  program  was  operated  with  a  mean  July  climate  for  California  that  speci¬ 
fied  P(h),  T(h),  and  v(h)  profiles  over  the  range  h  *  3.8  to  30  km.  Dry  air  zenith 
attenuation  Ap  and  a  water  vapor  attenuation  slops  ay  s  Ay/V  were  calculated  by  numer¬ 
ical  integration  and  listed  in  Table  7.  Values  of  Ap  have  been  converted  into  bright¬ 
ness  temperatures  Tp  by  assuming  an  average  temperature  To  *  260  K  consistent  with 
the  T(h)  profile.  Brightness  data  below  10  GHz  are  particularly  sensitive  to  the 
correct  value  of  y0  in  the  dry  air  continuum  (12a).  Above  10  GHz,  the  water  vapor 


slope  ay  amounts  to  meaningful  values. 


At  90  GHz,  the  data  potnt  (Tx  -  Tp)/33  = 


1.70  K/rnm  was  used  as  a  reference  water  vapor  slope  zy  to  define  an  Hequivalent" 
temperature  of  305  K  for  conversions  from  path  attenuation  to  medium  brightness  at 
other  frequencies.  This  way,  all  six  Tg-data  in  Table  7  are  apportioned  in  a  con¬ 
sistent  manner  to  dry  and  wet  attenuation  terms  based  upon  MPM  predictions. 

in  nummary,  good  agreement  between  predicted  and  reported  responses,  foremost 
specific  attenuation,  was  found  over  a  wide  range  of  parameters- -test  frequencies 
varied  between  2.5  and  430  GHz  and  meteorological  conditions  as  follows: 

P  =  70  to  101  kPa.  T  =  -1C® to  Jb^C,  v  =  0  to  25  g/a*  and  RM  =  0  to  99X.  All  in  all, 
some  credibility  can  be  given  to  the  computer  program  MPM  detailed  in  Section  2.1. 
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